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ABSTRACT 

We conducted a study of the environment around the supernova remnant CTB109. We found that 
the SNR is part of a large complex of H II regions extending over an area of 400 pc along the Galactic 
plane at a distance of about 3 kpc at the closer edge of the Perseus spiral arm. At this distance CTB109 
has a diameter of about 24 pc. We demonstrated that including spiral shocks in the distance estimation 
is an ultimate requirement to determine reliable distances to objects located in the Perseus arm. The 
most likely explanation for the high concentration of H II regions and SNRs is that the star formation 
in this part of the Perseus arm is triggered by the spiral shock. 

Subject headings: circumstellar matter, H II regions, ISM: clouds, ISM: kinematics and dynamics, 
supernova remnants 



1. INTRODUCTION 

The supernova remnant (SNR) CTB109 was discovered 
by Gregory & Fahlman (1980) at X-rays with the Einstein 
Satellite. It is believed that the peculiar semicircular shape 
of the SNR results from its interaction with a giant molecu- 
lar cloud (Tatematsu et al. 1987) which inhibits its expan- 
sion to the west. It has an unusual double-shell structure 
in radio, with an inner half-shell surrounded by a wider 
outer shell; both shells have similar radio surface bright- 
ness at 1420 MHz (see Fig. 1). CTB109 was classified as 
a shell-type SNR based on its radio properties (Downes 
1983). The supernova explosion left a pulsar behind that 
is observable in X-rays (Fahlman & Gregory 1981) but no 
radio emission has been detected from the pulsar and there 
is no evidence for a pulsar wind nebula. 

Establishing the distance to a SNR is usually quite dif- 
ficult, and CTB109 seems to offer a particular challenge. 
£ - D relations were used by Sofue et al. (1983) and Hughes 
et al. (1984) to derive distances of 4.1 kpc and 5.6 kpc re- 
spectively. The newest approach to £ - D is that of Case & 
Bhattacharya (1998), who use only shell-type SNRs; this 
gives a distance of 4.7 kpc. However, the S-D technique, 
while the oldest, is probably the most disputed method 
of distance determination for SNRs. There are many sta- 
tistical studies which show that there is no evidence for 
a working E-D relation (Green 1984; Berkhuijsen 1986, 
1987). 

The spectroscopic distance of H II regions close to 
CTB109 offers an alternative approach, and Tatematsu 
et al. (1990) quote values of 3.6 - 5.4 kpc . However, 
these distances are based on very old measurements and 
only two nearby H II regions were used in this publication. 
Despite the results mentioned above, all previous authors 
have used a distance of 4 kpc for CTB109. 

The primary motivation for the present study is to de- 
termine a more reliable distance to this SNR. We use spec- 
troscopic distances and radial velocities of 16 nearby H II 
regions and compare these values with those obtained for 



CTB109. We also include a comparison of absorption pro- 
files of the radio bright H II regions with foreground H I 
column densities of CTB109 calculated from X-ray absorp- 
tion and discuss the cold environment of the SNR based 
on high-resolution H I and CO measurements from the 
Canadian Galactic Plane Survey (CGPS). 

We have also obtained 408 MHz and 1420 MHz radio 
continuum data including polarization at 1420 MHz. A 
joint spectral index and polarization analysis of these data 
in conjunction with 2.7 GHz, 4.85 GHz, and 10.55 GHz 
measurements from the Effelsberg 100 m radio telescope 
will be published elsewhere (Kothes et al. 2002). 

2. OBSERVATIONS AND DATA ANALYSIS 

The CGPS is described in detail by Taylor et al. (2002). 
To study the environment of the SNR CTB109 we have 
used parts of the CGPS which derive from observations 
with the Synthesis Telescope of the Dominion Radio As- 
trophysical Observatory (DRAO) (Landecker et al. 2000) 
and CO observations (Heyer et al. 1998) from the Five 
College Radio Astronomy Observatory (FCRAO). For the 
DRAO data angular resolution varies as cosec (declination) 
and therefore changes slowly across the final maps. Pa- 
rameters for all observations used can be found in Tab. 
1. To ensure accurate representation to the largest scales 
for the DRAO data we have incorporated single-antenna 
observations after suitable filtering in the Fourier domain. 
Single-antenna H I data were obtained from a survey of 
the CGPS region made with the DRAO 26-m Telescope 
(Higgs & Tapping 2000). The single- antenna data for the 
1420 MHz radio continuum were obtained from the Effels- 
berg 1420 MHz Survey of the Galactic Plane (Reich et al. 
1997). 
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3. THE DISTANCE TO CTB109 

In light of all the distance problems around CTB109, we 
have embarked on a more thorough analysis and have tried 
to determine a more accurate distance to this SNR. Aban- 
doning the S-D relation, the most commonly used method 
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to determine the distance to a SNR is to deduce the ra- 
dial velocity of associated neutral hydrogen or molecular 
material. A Galactic rotation model is then used to as- 
sign to the radial velocity a distance from the Sun. The 
molecular material associated with CTB109 has a radial 
velocity of about —51 km/s (Tatematsu et al. 1990) (see 
also section 3.2), implying a kinematic distance of 5 kpc 
(using a flat rotation model with vq = 220 km/s and Rq = 
8.5 kpc). However, the validity of the fiat rotation model 
has been challenged for the Perseus arm. Roberts (1972) 
pointed out that H II regions show a high discrepancy 
between their kinematic and spectroscopic distances. Fur- 
thermore, O-type stars in the Perseus arm region often 
show interstellar absorption lines, arising of course in ma- 
terial in front of them, with velocities more negative than 
the radial velocity of the star itself. According to the flat 
rotation model this material should be behind those stars. 
Roberts (1972) successfully modeled these discrepancies 
with the two-armed spiral shock (TASS) model. 

In Tab. 2 we list parameters of 16 Sharpless H II re- 
gions between Galactic longitudes 104° and 114° taken 
from Brand & Blitz (1993, and references therein), an ex- 
tensive study of the velocity field of the outer Galaxy. In 
Fig. 2 the radial velocity of these H II regions is plot- 
ted as a function of their spectroscopic distances. There 
are two groups of H II regions, separated by their ve- 
locity and distance. The first group comprises sources 
with distances well below 2 kpc and radial velocities be- 
tween and —20 km/s. These sources are located in the 
Orion Spur, and their spectroscopic and kinematic dis- 
tances agree quite well. There is a large gap which corre- 
sponds to the interarm region between the Orion Spur and 
the Perseus arm. The H II regions of the second group are 
all located in the Perseus arm, and most are much closer 
than predicted by the flat rotation curve. This result fits 
well with Roberts' predictions (Roberts 1972, Fig. 4). Ac- 
cording to Roberts' model, the radial velocity in the di- 
rection of CTB109 would drop to about —50 to —55 km/s 
at the position of the spiral shock at a distance of about 
2.5 kpc. It then slowly rises until it rejoins the flat rota- 
tion curve at about 3.5 kpc. The only H II regions with 
velocities between —40 and —60 km/s whose spectroscopic 
distances exceed 4 kpc are the Shl47/8/9 complex and 
Shi 56 (and these two distances have the biggest errors). 
Independent distance estimates for these two H II regions 
obtained from infrared measurements by Wouterloot et al. 
(1988) are about 3.5 kpc for both. The distance ambiguity 
for Perseus arm velocities is also indicated by the presence 
of many H I self absorption features (HISA) , found by Gib- 
son et al. (2000). The HISA phenomenon requires warm 
neutral gas behind cold absorbing gas. The emission of 
the background hydrogen is then absorbed by colder fore- 
ground material at the same radial velocity. This can only 
occur when two distances correspond to the same radial 
velocity: a flat rotation curve does not permit this in the 
outer Galaxy (e.g. Fig. 2), but the Roberts TASS model 
does. 

We have also produced absorption profiles for the radio- 
bright H II regions in our list. From these we were able to 
calculate accurate foreground H I column densities, also 
listed in Tab. 2. These values are comparable to those ob- 
tained by Rho & Petre (1997) towards CTB109 from X-ray 



absorption. They found absorbing H I column densities 
between 8 and 10 • 10 21 cm" 2 . Patel et al. (2001) derived 
an absorbing H I column density of N# = 9.3 • 10 21 cm~ 2 
for the central pulsar from Chandra observations. For 
the radio-bright southern spot which is absorbed by fore- 
ground material in our data we measure an absorbing H I 
column density of 8.4 ± 0.8 • 10 21 cm~ 2 . 

The absorption profiles obtained for the radio bright 
southern spot of CTB109, the two closest H II regions 
Shl49 and Shl52, and a nearby extragalactic source are 
plotted in Fig. 3. All profiles are very similar for the local 
gas. The profiles of the three Galactic objects are also very 
similar for the Perseus arm gas, indicating that they are 
located at comparable distances. However, the extragalac- 
tic source has an additional component around —45 km/s 
which is completely missing in all of the Galactic sources 
(and so must arise beyond the Galactic sources). This is 
another indication for the existence of the TASS, because 
the radial velocity within the Perseus arm would be in- 
creasing with distance in the TASS model while it would 
be decreasing if the flat rotation curve applied. This evi- 
dence also favors the closer distance for Sh 149, since ma- 
terial with a radial velocity of around -45 km/s, which is 
absorbed by the extragalactic source but not by Sh 149, 
should be located at around 4 kpc according to Roberts 
(1972). We should note at this point that the deduction 
of the systemic velocity of the Galactic sources from their 
absorption spectra is not possible due to the peculiar be- 
havior of the rotation curve in the Perseus arm. 

CTB109 contains a pulsar, discovered by Fahlman & 
Gregory (1981), and evidently is the product of the explo- 
sion of a massive star. Hence it seems quite reasonable 
to associate it with the other massive stars in the vicin- 
ity, which are exciting the nearby H II regions. In light of 
all this new information we propose that CTB109 and the 
nearby H II regions are part of the same complex, which 
is located in the shock zone of the Perseus arm as de- 
scribed by Roberts (1972). This would imply a distance of 
3.0 ± 0.5 kpc. The mean distance of all H II regions listed 
in Tab. 2, weighted by 1/cr 2 , is 3.1 ± 0.2 kpc, which sup- 
ports the above mentioned distance estimate. The linear 
size of the SNR at 3.0 kpc would be 24 pc. The existence 
of this group of H II regions at about the same distance 
suggests that star formation in the Perseus arm in this 
vicinity was triggered by the spiral shock. 

4. THE COLD ENVIRONMENT OF CTB109 

The 1420 MHz radio continuum image of CTB109 taken 
from the CGPS is shown in Fig. 1. The radio emis- 
sion shows two shells of radius 8' and 18'. Both shells 
are incomplete in the northwest where the radio emission 
abruptly stops. The binary pulsar IE 2259+586 lies near 
the geometric center of the larger shell. There is another 
bright emission structure in the south where both shells 
overlap. The remnant is brightest to the northeast and 
southeast with an intervening depression. A striking fea- 
ture of the radio emission is the presence of several thin 
filaments between the shells and even some on top of them, 
indicating complex dynamical structures within the rem- 
nant and around it. 

4.1. The molecular material 
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The molecular environment of CTB109 was studied in- 
tensively by Tatematsu et al. (1987, 1990) on the basis 
of CO observations with the Nagoya 4-m radio telescope 
and the 45-m telescope at the Nobeyama Radio Observa- 
tory, respectively. They described the massive molecular 
cloud to the west of CTB109 and found a ridge of molec- 
ular material extending inside the SNR, showing an anti- 
correlation with the X-ray emission. They found no evi- 
dence of shocked CO inside or around the remnant. In Fig. 
4 we plot channel maps of our CO data towards CTB109 
around the systemic velocity of the SNR. We should note 
at this point that, in contrast to the customary Perseus 
arm velocity profile, in this vicinity higher negative veloc- 
ity indicates a shorter distance (see also section 3). 

The massive molecular cloud to the west of the SNR is 
very prominent in the channel maps. This cloud blocks 
the shock wave from further expansion in this direction. 
According to Tatematsu et al. (1987) this cloud is respon- 
sible for the peculiar shape of the SNR. In the northern 
part there is an extension from the cloud into the SNR 
separating both radio shells. This CO ridge appears at 
about —46 km/s to —47 km/s and vanishes at —52 km/s 
while the main body of the cloud appears at a slightly 
more negative velocity, moving in gradually from the west 
until it disappears in the same direction at —54 km/s to 
—55 km/s. This indicates that the main body of the cloud 
is located closer to us than the ridge, or that the CO ridge 
is moving away from us relative to the big cloud. It also 
shows that the core of the massive cloud is to the west 
and the SNR shock wave has been reflected by its eastern 
edge. The dark cloud at the tip of the ridge, which is more 
pronounced in the integrated map in Fig. 1 , causes strong 
absorption in the X-ray data (Tatematsu et al. 1987, 1990; 
Rho & Petre 1997) indicating that it lies in front of the 
SNR or is associated with it. The only possible radial 
velocity for CTB109 would be between —50 km/s and 
—52 km/s where both structures, the ridge and the mas- 
sive cloud, are present. This also explains the weak outer 
radio shell beyond the CO ridge. The shock wave traveling 
in that direction has mostly been blocked by this ridge. 

4.2. Neutral hydrogen 

Associated neutral hydrogen features are less pro- 
nounced than those seen in CO, but the general structure 
seen in the CO data is also present. The massive cloud 
to the west extends further to the south and is generally 
closer to the remnant. The ridge in the north is also vis- 
ible, although it is seen slightly further south, closer to 
the inner SNR shell, and can be followed further around 
the inner shell. It appears to be smoother, which makes 
it difficult to detect in the channel maps (Fig. 5), but 
in the integrated map of Fig. 1 it is quite obvious. The 
SNR seems to be located at a density gradient in the H I 
distribution most prominent in the three channel maps in 
the second row of Fig. 5. This gradient goes from high 
densities in the west smoothly down to a hole in the H I 
distribution cast of CTB109. Fig. 1 also reveals a diffuse 
shell of HI surrounding almost the whole outer radio shell 
of the SNR. Note that the deep depression in the area be- 
tween the bright spot in the south and the northeastern 
part of the outer shell (in the H I channel maps between 
-47 and -53 km/s) cannot be caused by absorption since 



the continuum emission of the outer shell in this area is 
not strong enough. Therefore this relative absence of H I 
emission is due to the presence of the smooth outer shell. 
Absorption of the bright southern spot is also detectable 
in the channel maps. 

5. DISCUSSION 
5.1. CTB109 as part of a large complex 

The distribution of H II regions and SNRs in the vicinity 
of CTB109 is displayed in Fig. 6 in relation to the Perseus 
arm locations from Roberts (1972) and Taylor & Cordes 
(1993). Apparently the kinematic distance for CTB109 
derived from the flat rotation curve would place the SNR 
in the interarm region behind the Perseus arm. Even the 
4 kpc distance assumed by various earlier authors, would 
be outside the spiral arm. Since the SNR is associated with 
dense molecular and atomic material an interarm location 
would be very unlikely. 

In the area shown in Fig. 6 there are two concentrations 
of sources within the Perseus arm. The first is around 
CTB109 and the supernova remnant Cas A and the sec- 
ond is around the Tycho SNR. There is a third concen- 
tration at higher longitudes around the W3/4/5 complex 
and the supernova remnant HB3 (not shown here). Ob- 
viously, while traveling through the Perseus arm the spi- 
ral shock triggered star formation in high-density regions, 
where we now find these concentrations of young stellar 
objects, indicated by the compact H II regions and su- 
pernova remnants which are the result of the explosion of 
massive young stars as well. 

5.2. Neutral material in the vicinity of CTB109 

CTB109 is located at a density gradient which goes from 
very high densities on the molecular cloud in the west to a 
gap in the emission in the east of the SNR. Apparently 
the progenitor star was formed at the edge of a dense 
molecular cloud while the star in Sh 152 was formed in 
its centre. To the west, the shock wave expanded into 
the dense cloud which decelerated it very quickly. To the 
east, however, the shock wave is expanding into a moder- 
ately dense medium. The structure of the surrounding H I 
indicates that the SNR is not expanding inside a stellar 
wind bubble, because we would expect such a structure to 
have a more pronounced outer boundary. It is more likely 
that the missing emission from the interior is the result 
of taking away the neutral hydrogen by ionizing it with 
the expanding shock wave. This would imply a progenitor 
star of type B2/3 since more massive stars would have a 
strong stellar wind. 

Between the outer and inner shell in the north we find 
the ridge of cold material separating both structures. The 
fact that the H I is smoother and mostly concentrated 
outside the thin, dense molecular ridge indicates an inter- 
action with the SNR in which the surface of the molec- 
ular structure was dissociated or even evaporated by the 
expanding shock wave. The lack of X-ray emission coin- 
ciding with the ridge indicates that the ridge is located on 
the near side of the remnant and is absorbing the X-ray 
emission which originates behind it. 

The outer radio shell shows two very prominent features. 
These are the bright knot to the south and the bright part 
of the northeastern shell. Both features are located within 
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bright parts of the surrounding H I (see Fig. 1) indicating 
that these parts of the remnant are expanding into a higher 
density medium than the part between them where the 
emission of the radio shell is weakest and the surrounding 
H I seems to have almost a gap. This suggests a close rela- 
tion between the radio brightness of the SNR's expanding 
shell and the density of the medium it is expanding into, 
which is of course expected. 

6. SUMMARY 

We have presented new radio continuum data for the 
SNR CTB109 together with H I and CO observations of 
the surrounding medium. Our data show that the radio 
continuum emission from the remnant is closely related to 
the surrounding cold material. Analysis of the H I and CO 
dynamics and comparison of the results with parameters 
of nearby H II regions has led us to the conclusion that the 
SNR is located at a distance of 3.0 ± 0.5 kpc, as opposed 
to larger distances previously published in the literature. 
It also implies that the SNR, together with several H II 



regions, is part of a large complex created by the spiral 
shock present in the Perseus arm. We have shown that 
inclusion of the spiral shock is necessary to produce reli- 
able kinematic distances in this part of the Galaxy, and 
we have demonstrated that, in seeking the distance to a 
supernova remnant, it is important to take a broader view 
of the environment. The wide-field, high-resolution data 
of the CGPS, revealing multiple components of the ISM, 
are ideal for a study of this type. 

We like to thank Tom Landeckcr for careful reading of 
the manuscript. We also like to thank Christopher Brunt 
for many useful discussions about the spiral shock in the 
Perseus arm. The Dominion Radio Astrophysical Obser- 
vatory is a National Facility operated by the National Re- 
search Council. The Canadian Galactic Plane Survey is a 
Canadian project with international partners, and is sup- 
ported by the Natural Sciences and Engineering Research 
Council (NSERC). 
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Table 1 



CGPS Data Parameters at the position of CTB109 


Parameter 


Value 


Pixel Size 


18" 


Number of Spectral Channels 


256 


Channel Separation 


0.8243 km s" 1 


LSR Velocity Range 


+45 km/s to -170 km/s 


H I Spectral Resolution 


1.319 km/s 


CO Spectral Resolution 


1.0 km/s 


H I Beam Size a 


68" x 59" at -63° 


CO Beam Size 


45" x 45" 


H I RMS Noise per Channel 


2.5 K 


CO RMS Noise per Channel 


0.2 K 



a Position angle of the beam is with respect to the Galactic 
pole. 



Table 2 

Parameters of 16 nearby H II regions taken from Brand & Blitz (1993). The errors in the H I column density 

VALUES ARE BETWEEN 10 AND 20 %. 



Name Gal.-Coor. Distance [kpc] vlsr [km/s] N H i [10 21 cm 2 ] 



Shl35 


104 


.59+1.37 


1.4+0.4 


-20 


.7+0.5 




Shl37 


105 


.15+7.12 


0.6+0.2 


-10 


.3+1.4 




Shl39 


105 


.77-0.15 


3.3+1.1 


-46 


.5+0.5 




Shl40 


106 


.81+5.31 


0.9+0.1 


-8 


.5+1.0 




Shl42 


107 


.28-0.90 


3.4+0.3 


-41 


.0+0.5 




Shl49 


108 


.34-1.12 


5.4+1.7 


-53 


.1+1.3 


8.4 


Shl52 


108 


.75-0.93 


3.6+1.1 


-50 


.4+0.5 


8.2 


Shl54 


109 


.17+1.47 


1.4+0.4 


-11 


.5+0.9 




Shl55 


110 


.22+2.55 


0.7+0.1 


-10 


.0+1.5 




Shl56 


110 


.11+0.05 


6.4+2.0 


-51 


.0+2.0 


8.7 


Shl57 


111 


.28-0.66 


2.5+0.4 


-43 


.0+2.0 




Shl58 


111 


.54+0.78 


2.8+0.9 


-56 


.1+1.1 




Shl59 


111 


.61+0.37 


3.1+1.2 


-56 


.0+1.0 


11.0 


Shl61B 


111 


.89+0.88 


2.8+0.9 


-51 


.9+0.7 




Shl62 


112 


.19+0.22 


3.5+1.1 


-44 


.7+0.5 


7.7 


Shl63 


113 


.52-0.57 


2.3+0.7 


-44 


.9+3.8 
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Fig. 1. — upper left: Radio continuum image of the SNR CTB109 taken from the CGPS at 1420 MHz. The shading goes from 5 K (white) 
to 18 K (black). The position of the X-ray pulsar 1E2259+586 is marked by a white asterisk. Upper right: CO around CTB109 integrated 
between -45 km/s and -55 km/s. Contours are from 0.2 K (white) to 4.2 K (black) in steps of 0.5 K. Lower left: H I around CTB109 integrated 
between -45 km/s and -55 km/s. Contours arc from 70 K (black) to 94 K (white) in steps of 3 K. 



Fig. 2. — Radial velocities of H II regions in the vicinity of CTB109 as a function of their spectroscopic distances. The solid line represents 
the expected behaviour according to the flat rotation model with vq = 220 km/s and Rq = 8.5 kpc. The radial velocity of CTB109 is 
indicated by the dashed line. 



Fig. 3. — Absorption and emission profiles of the southern bright spot of CTB109, Shl49, Shl52, and the extragalactic radio source KR 68a. 



Fig. 4. — CO channel maps towards the SNR CTB109 taken from the CGPS. White contours indicate the radio continuum emission at 
1420 MHz. 



Fig. 5. — H I channel maps towards the SNR CTB109 taken from the CGPS. White contours indicate the radio continuum emission at 
1420 MHz. 



Fig. 6. — H II regions (open circles) and SNRs (open stars) in the vicinity of CTB109. The Perseus arm as defined by Roberts (1972) is 
indicated by the dashed lines. The location of the maximum electron density in the Perseus arm as determined by Taylor & Cordes (1993) is 
marked by the dotted line. The names of the individual SNRs are noted. For Shl56 and the Shl47/8/9 complex the distances from Woutcrloot 
et al. (1988) are used and for all other H II regions the spectroscopic distances listed in Brand & Blitz (1993). The distances for the SNRs 
Cas A and Tycho were taken from Reed et al. (1995) and Chevalier ct al. (1980), respectively. 
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